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We derive an exact expression for entropy production during effusion of an ideal gas driven by 
momentum transfer in addition to energy and particle flux. Following the treatment in Phys. Rev. 
E 74, 021117 (2006), we construct a master equation formulation of the process and explicitly verify 
the thermodynamic fluctuation theorem, thereby directly exhibiting its extended applicability to 
particle flows and hence to hydrodynamic systems. 

PACS numbers: 05.70.Ln, 05.40.-a, 05.20.-y 



43 

9 L,: 

CZ2 
Jh ■ 

Td ; 

CZ3 ■ 

CD ■ 
• i— i . 

43 ' 



> 

(N 

O 

O 

>>: 

43 

Or 



X 



I. INTRODUCTION 

Since the pioneering work of Onsager fl[ on the relation 
between linear response and equilibrium fluctuations, his 
insights have been further formalized in, for example, 
the theory of linear irreversible processes Q and the 
fluctuation-dissipation theorem [H . Over the past decade 
some new surprising results have been discovered that 
suggest relations valid far awayfrom equilibrium, notably 
the fluctuation |4( and work [5j theorems. The fluctua- 
tion theorem, originally demonstrated for noncquilibrium 
steady states in thermostated systems, has been proven 
in a number of different settings. Basically, it states that 
during an experiment of duration t, it is exponentially 
more likely to observe a positive entropy production AS 
rather than an equally large negative one, 



Pt(AS) 
Pt(-AS) 



= e AS/k . 



(1) 



In the application to nonequilibrium steady states, the 
above result is typically only valid in the asymptotic limit 
t — > oo and expresses a symmetry property of large devi- 
ations. 

We address another scenario in which the system is 
perturbed out of a state which is initially at equilibrium. 
The so-called transient fluctuation theorem is then valid 
for all times t. We consider the problem of a Knudscn 
flow between ideal gases that have overall non-zero mo- 
mentum. In this case, the stationary state is reached in- 
stantaneously, so that there is no distinction between the 
transient and steady state versions of the theorem. We 
show that the system obeys a detailed fluctuation theo- 
rem which includes |T]) as a special case. Our calculation 
is an extension of the one given in Q to include momen- 
tum transfer. The interest of this extension is manifold. 
First, the derivation of the fluctuation theorem is some- 
what more complicated since the momentum is a quantity 
which is odd under velocity inversion. Second, momen- 
tum, together with particle number and energy, are the 



conserved quantities whose transport forms the basis of 
hydrodynamics. Our derivation therefore puts the fluctu- 
ation theorem fully in this context (see also Q ) • Finally, 
as a bi-product, we calculate the Onsager matrix for the 
Knudsen flow problem including momentum transport. 

In Sec. |TT] we formulate the fluctuation theorem for effu- 
sion with momentum transfer for an ideal gas. Section lTTTl 
generalizes the derivation of the master equation and cu- 
mulant generating function in Q to the case with mo- 
mentum transfer. Verification of the fluctuation theorem 
is detailed in Sec. lIVi and the lowest order cumulants are 
exhibited in Sec|V] We use these results to verify the On- 
sager relations for this nonequilibrium system in Sec I VII 
We end with a brief conclusion in Sec. IVIIl Some details 
of the calculations are presented in appendices. 



II. FLUCTUATION THEOREM FOR EFFUSION 
WITH MOMENTUM TRANSFER 

We begin by considering two (infinitely) large neigh- 
boring reservoirs, A and B, each of which contains an 
ideal gas of uniform density in equilibrium at temper- 
ature Ti, i € A, B. In addition, the particles of gas i have 
an overall center of mass velocity Vi in the e x direction 
(Fig. [T|). That is, the velocity distributions of the gas 
particles take the Maxwellian form 



2nkTi 



3/2 



exp 



m(v - Vje x ) u 
2kf t 



(2) 



The two reservoirs are separated by a common adiabatic 
wall parallel to the e x direction, with a hole of surface 
area a whose linear dimensions are small compared with 
the mean free path of the particles. As a result, the 
local equilibrium in each reservoir is not disturbed by the 
exchange of mass, heat, and momentum during a finite 
time interval t in which the hole is open. Upon a transfer 
of total energy At/, particles A7V, and momentum Ap x 
during this time interval, the overall change in entropy 
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FIG. 1: The ideal gas in compartment A is characterized by 
temperaturerTA , number density ha, and flow velocity Va- 
Similarly, the ideal gas in compartment B is characterized by 
temperature Tb, number density jib, and flow velocity Vb- 
Velocities Va and Vb are represented by the vertical arrows 
within each compartment and are taken to be parallel to the 
adiabatic wall separating the the two gases. 



for the system is given by 
AS = AS A + AS B 



-±AU- 
T A 


Ta 


mV\ 
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)AN- 


\- —Ap x 
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Tb 


mVl 
2T B 
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- 7^A Px 
-Lb 



(3) 



= AuAU + A N AN + A Px A Px , 

where we have introduced the thermodynamic forces for 
energy, particle, and momentum transfer 



1 

f B 
T A 



An = =- - =- 



A — ^ A 

An — 



= /clog 



Va 
T a 



1 

T~a 
rnV~A 
~Wa~ 

riA ( T B 
n B 
Vb 
T b 



_ /Ms _ raVg 
Tb 2T b 



T a 



I ( mVj 
[ 2T R 



~ 7 2T A ~' 



(4) 



In the equation for An, we have used the expression for 
the chemical potential ll of an ideal gas at rest. Since the 
explicit expression for the thermodynamic forces in sys- 
tems with momentum is not readily available, we provide 
a brief derivation in Appendix El 

The variables AU, AN, and Ap x all correspond to fluc- 
tuating quantities influenced by single particle crossings 
on each side of the adiabatic wall. As a result, the total 
entropy production AS*, which will be observed during 



the time duration t, is likewise a fluctuating quantity. 
However, time reversal symmetry of the microscopic dy- 
namics a relation for the probability distribution of this 
entropy production, as expressed in Eq. ([1]). Due to the 
absence of memory effects, AS is in fact a stochastic pro- 
cess with independent increments: contributions to AS 
from any two equal, non-overlapping time intervals are 
independent identically distributed random variables. It 
is therefore convenient to introduce the cumulant gener- 
ating function, which takes the form 



-tfl(A) 



(5) 



The fluctuation theorem, Eq. [T] implies the following 
symmetry property: 



g(X) =g(k- 1 - A). 



(6) 



As the derivation in Appendix [B] points out, one can 
augment the observation of the entropy production with 
additional variables, while retaining the form of the fluc- 
tuation theorem. Hence the following more detailed fluc- 
tuation theorem, which is expressed in terms of the joint 
probability density involving all three conserved quanti- 
ties, particle number, momentum and energy, is obtained: 



P t {AU,AN,A Px ) 
P t (-AU,-AN,-A Px ) 



= e AS/k . 



(7) 



Since the increments of AU, AN, and Ap x are also inde- 
pendent, we can write the corresponding cumulant gen- 
erating function as 



-{XuAU+\ N AN+\ Px A P:c )^ _ e -tg(\u ,X N ,A P J 



(8) 



The detailed fluctuation theorem then requires the fol- 
lowing symmetry relation, similar to Eq. ([B]): 

<?(A[/,Ajv, A Px ) 

= g{A u k^ 1 - \u,Ank~ l - AAr,„4p x fc -1 - A P J. 

(9) 

Note finally that Eq. |J7J| , apart from implying the nor- 
mal fluctuation theorem Eq. |T]), also implies fluctua- 
tion theorems for particle, energy, and momentum trans- 
fer individually when the complementary thermodynamic 
forces are zero: 



Vt(AU) 
V t (-AU) 

Vt(AN) 
Vt(-AN) 

Vt{A Px ) 
V t {-Ap x ) 



e AS / fe ; An = A 
e AS / k ; Au=A 



, -0, 
= 0, 



(10) 



= e AS / fe ; An = A N = 



III. 



MASTER EQUATION AND CUMULANT 
GENERATING FUNCTION 



If we choose a sufficiently small time interval dt, the 
contributions to the quantities AU, AN, and Ap x arise 
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from individual particles crossing the hole. The kinetic 
theory of gases allows us to calculate the probability per 
unit time, Ti->j(E,p x ), to observe a particle with kinetic 
energy E = ^mv 2 and momentum p x = mv x crossing 
the hole from reservoir i to reservoir j. Specifically, the 
transition rate in question is given by (see Appendix [A"|) : 



an. 



x exp 



m(fc 7 rT i ) 3 / 2 
rn 



E - 



2 \ 1/2 

Px 

2m 



2kTj \ m m 



(11) 

with = {A, B) or — (B,A). Hence, the prob- 
ability density Pt(AU, AN, Ap x ) obeys the master equa- 
tion 

d t P t (AU,AN,Ap x ) 



dp x / dE T A ->b 

) Jp^/2m 

xP t (AU-E,AN-l,A Px - Px 



+ 



dp x / dE T b ^a 



(12) 



x P t (AU + E,AN+l,A Px + p x ) 

/oo poo 
dp x / dE 
-oo Jp^/2m 

x (T A ^ B + T b ^a) • 



We have written Ti—,j(E,p x ) without the arguments for 
economy of notation. We can take advantage of the con- 
volution structure of the integral operators by consider- 
ing the equation in Fourier space; that is, we multiply 
both sides of the equation by exp[— (XjjAU + AjyAiV + 
Xp^Apx)] and integrate AU, Ap x over all space and sum 
over all integers AN. We arrive at the expression 



d t P(Xu, Ajv,A p J = P(\u,\ N , ApJ/i 



+ P(^u, An, A Px )/ 2 - P(Xu,Xn, \*)h, 



(13) 



where 
h=e 

I 2 = e 

h = 



d Px dET A ^ B e~ {XuE+x ^\ 

oo Jp2/2m 
poo 

dp x dET B ^ A e {XuE+x ^\ 

' ''Px/ 2 ™ 



dp x 



p|/2m 



dE (E 



A^B 



T, 



B^A 



From this expression we can write g{Xjj, Ajv, A Px 
fined in Eq. ([8]), as 

g(X v , X N , A Px ) =I 3 - (h +h)- 



(14) 
, de- 



(15) 



The integrals (|14|1 can easily be performed by switching 
to the variable z = E—p 2 /2m and integrating z from zero 



to infinity and p x over all space, as before. We thereby 
arrive at our final expression for g(Xu, Xn, A Px ): 



](Xjj,Xn, A Px ) — a 



2-Km 



uaTT ( 1 



n B T B /2 ( 1 



1/2 



Ga 



(1 + kT A Xuf 
Gb 



(16) 



where 



G A = exp -X N 



Gb = exp A 



JV 



(l-kT B Xu)\ 

VlXu - kmT A Xl x + 2mV A X Vx 
2(1 + kT A Xu) 

mV^Xu + kmT B X 2 v ^ + 2mV B X Pjs \ 
2(1 - kT B X v ) I 



(17) 

Notice that g(Xu, Xn, A Px ) can be written as a sum of 
two contributions, 

g(Xu, Aat, A Px ) = 3a(A;7, Xn, A Px ) + gsiXu, Xn, A Px ), 

(18) 

with 



gA(Xu, Xn, A Px ) = a 



2irm 



1/2 



x n A T A /2 [ 1 



G A 



(1 + kT A XuY 



ffs(A[/, Aat, A Px ) — g 



2irm 



1/2 



(19) 



x n B T B /2 ( 1 



Gb 

(1 - kT B Xu) 2 



This additivity property arises from the statistical inde- 
pendence of the fluxes from A — > B and B — > A. 



IV. FLUCTUATION SYMMETRY 

We now proceed to explicitly verify the symmetry re- 
lation @ and hence the fluctuation theorem. Concep- 
tually, one may understand this symmetry relation as 
follows. Under the symmetry operation T - that is, un- 
der the transformation given by the r.h.s. of Eq. © - 
the term containing the exponential in gA (which we call 
gA,i) becomes the corresponding term from gs (which 
we call gB,i) and similarly, the original gB,i term be- 
comes gA, l, thereby preserving the overall structure of g. 
Mathematically, we can express this as: 



= 9B,1, 

T[9ba] = 9a,i, 



(20) 



where 

g A ,i = <r (k2nm) 1/2 n A T A l2 {\ + kT A Xu)- 2 G A , (21) 
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and 



9b, i = o {k2nm) 1/2 n B T B /2 (l - kT B X v ) 2 Gi 



(22) 



Explicitly, we have 



l~[g A ,i] = & 



2nm 



1/2 



n A T A /2 (1 + kT A AuY 



x exp < log 



n A ( T B 



n B \T A 



3/2' 



A 



AT 



x exp 



x exp 



-mVjAu + kmT A Al x 
2(1 + kT A A v ) 

-2toVaA Px 



2(1 + kT A Au) 



2nm 



1/2 



x exp 



x exp 



\n B T A 



■> \ I V B m _ V A m 

- A 1 2feT B wrx 



mVjAu 



2A 



kmT A A^ - 2mV A \ 
2A 



where 



(23) 



Au = 



A„„ = 



1 



1 



kT B 
V A 



kT A 
V B 



kT A kT E 



A N = 



A = 



V 2 m 

— A 



V^m 



2kT B 2kT A 
^(l-kT B Xu). 



(24) 



N ■ 



Following simplification, this reduces to .95,1. A simi- 
lar result holds for T[g By i], and therefore the fluctuation 
theorem symmetry is satisfied. 



V. CUMULANTS 

The joint cumulant Kijk of power i in energy flux, j in 
particle flux, and k in momentum flux appears as a coeffi- 
cient in the Taylor expansion of the cumulant generating 
function, namely, 



.9,1 



( Ay , \n , \x ) = - - 



(-ly+^x^x^ 



i\j\k\ 



(25) 

While our expression for g A (Xjj, A at, X Pm ) allows us to cal- 
culate joint cumulants of any order, we here mention only 
the first order results, which are relevant for verifying the 



Onsager relations in the subsequent section: 
kioo - (AU) 



to 



1 \ V2 

k \ I 1/2 

A 



n A T 1 A /2 (2kT A 



2nm J V 2 
-n B T B /2 (2kT B 



koio = (AJV> 



(26) 



to 



2irm 
K001 = (&Px) 

km 



1/2 



n A T A < - n B T B ' , , 



to ' (n A V A T A /2 - n B V B T B /2 



Note that the cumulant associated with energy (AU) 
contains terms corresponding to both particle transport 
and momentum transport. 



VI. ONSAGER RELATIONS 

Averaging Eq. ([3]) and taking the time derivative leads 
us to an equation for the average entropy production, 



dt 



(AS) = JuAu + J N A N + J Px A Vx , (27) 



with the macroscopic fluxes Jx defined as 



~ n B T B < (2kT B + 

= *(^) 1/ \n A V A TT-n B V B Tr 

(28) 

While these fluxes are in general complicated nonlinear 
functions of the affinities (Au ,An ,A Pn .) , near equilibrium 
we can write: 

AT AT 
T A = T-— , T B = T+—, 

An An , nn . 

n A = n — , n B = n+— — , (M) 



2 ' 

Ay 

2 ' 



v A = v-=-, V B = V + 



2 
AV 

~2~' 



and expand the forces and fluxes to first order in the 
small deviations AT, An, and AV. To linear order, the 
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thermodynamic forces become where J = (Ju,Jn,J Pi ) T and A — (Air, An, A Px ) T . 

The Onsager matrix O is given by: 

AT 

An - - j^-; 

mV ,3k mV 2 . . _ fc . / oa n 

^ = — Ay+(--^)AT--An (30) 

•>i ~~ T T 2 ' 

Taylor expansions of the fluxes J7i/,j7iv,and [Eq. ([28])] 
allow us to write 

J = OA, (31) 
I 



2mn 



1/2 



nT 3 / 2 



6fcT + 3mV 2 



(mV 2 ) 2 
4kT 



2 + 



2kT 



i 

kT 



mV 
kT 



raV 
kT 

(mV) 2 
r kT 



(32) 



which clearly has the required symmetry Oij = Oji. The 
Onsager relations [Eq. (|3 lj) ] fully detail the complex cou- 
pling between energy particle, and momentum transport 
in the linear regime. 

Note that in the case of moving gases, V ^ 0, the 
presence of a temperature gradient alone (An = AV — 0) 
is sufficient to produce a nonzero net flux of momentum. 
Note also that when there is only a momentum gradient, 
AT = Art = 0, the heat, particle, and momentum fluxes 
reduce to: 

ak^nT^ , __. A Tr 

Ju = -— -r-(mV)AV, 

(2-Km) 2 

Jn = 0, (33) 

akinT? 
Jv = - Tr) -rmAV. 

Therefore, when we choose the velocities to be equal but 
opposite so that 1^ = 0, the only nonzero flux is due to 
momentum transport. In other words, momentum ex- 
change takes place without a net exchange of particles or 
energy. 

VII. CONCLUSION 

The work and fluctuation theorems are quite remark- 
able. They are basically one further step in Onsager's 
program to take into account the time-reversal symmetry 
of the microscopic dynamics. This results in a stringent 
constraint on the probability density of the entropy pro- 
duction. The implications of this result are still being 



explored. In this paper we have shown by an explicit mi- 
croscopically exact calculation that the fluctuation the- 
orem applies for the effusion between ideal gases with 
non-zero overall momentum. This sets the stage for the 
application of the formalism in fluctuating hydrodynam- 
ics. 
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APPENDIX A 

We can derive the thermodynamic forces for an ideal 
gas of N particles in volume V with non-zero momentum 
by considering the entropy S(U, N, V) of a gas at rest as 
a function of U, the total energy, N, and V. Because 
adding an overall velocity to the gas does not change its 
volume in phase space and hence its entropy we can write 
the entropy S(U,N,V,p) of a flowing gas which depends 
on momentum p in terms of the entropy Sq(U, N, V) of a 
gas at rest: 

S(U, N, V,p) = S ((U - ^—), N, V) - S (e, N, V), 

( A1 ) 

where e = U — p 2 /2Nm represents the internal energy 
of the gas. The Sackur- Tetrode formula Q provides the 
explicit expression for S(e,N,V) which, with Eq. (|A1[) 
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leads to 



S = kN log 



such a realization is then given by 




3/2' 



2mN 



N 



-kN 



log 



\ 3h 2 
(A2) 



Here h is Planck's constant and m is the mass of a single 
gas particle. We can write the total entropy change of 
the effusion process considered here as 



dS 



dS A 
dU A 
dS B 



dU A 



dU{ 



dp( 



dU B 

dS B 
dU B 
dS R 



dS A 
dN A 
dS B 
dN B 

dU A 



dN A 



dN B 



dN( 



dS A 
dp A 
dS B 
dp B 

dS B 

dN R dN A 



dp A 
dp B 



dS / 



(A3) 



OSa 

dp B dp A 



), 



where Si corresponds to Eq. (|A2j) with U — ► Ui, N — > Ni, 
and p —>■ Pi, i E A, B, and we have used momentum, 
energy, and particle conservation to write dU — —dU A = 
dU B , dN = —dN A — dN B , and dp = —dp A = dp B . 
Performing the above calculations and considering that 
the total energy U of an ideal gas with overall momentum 
p at temperature T is given by 



U = 



3N p 2 
2kf + 2m' 



(A4) 



we arrive after simplification at the expressions given in 
Eq. 0. 



APPENDIX B 

We give a derivation of the fluctuation theorem, 
Eq. ([7]), by adapting to the present case the procedure 
introduced in [lOj. We consider the Hamiltonian evolu- 
tion of a system, consisting of two disjoint subsystems 
A and B initially at equilibrium characterized by micro- 
canonical distributions with total particle number, mo- 
mentum, energy and volume equal to Ni, pi, Ui and Vi, 
i = A, B, respectively. At the initial time, the constraint 
separating both systems is broken. It is assumed that 
this can be achieved without any external work, mo- 
mentum exchange or other perturbation of the subsys- 
tems. This is clearly the case for the opening of a hole in 
the adiabatic wall separating ideal g clSGS, cLS considered 
here. After a time interval of duration t, the constraint 
is again introduced at no cost of energy or momentum. 
One records the new values of the parameters N-, p[ and 
U[. The amounts (AN, Ap, AU) that are transported 
from system A to system B will depend on the specific 
run, i.e., on the starting configuration at t = 0. Let the 
volume in phase space corresponding to the initial states 
that lead to the transport of these amounts be denoted 
by f2(jVi,pi ,Ui) (AN, Ap, AU). The probability to observe 



{Ni, Pi ,Ui) 



.)(AN,Ap,AU) 



9. 



P (NitPitUi) (AN,Ap,AU) 

(Bl) 

where ^(Ar^pi.c/i) is the total phase space volume. We 
now apply this very same result for parameter values 
N-, —p'i and U[, i = A,B, and consider the probability 
of transporting the amounts (— AN, Ap, — AU). Apart 
from velocity inversion, the final values in this proce- 
dure are then the initial ones of the first scenario, i.e., 
(Ni, —pi, Ui). The corresponding probability reads 

P (N >,- p >.,u<)(-AN,Ap,-AU) 

= n {Nl ,- p[ , UD (-AN,Ap,-AU) (B2) 

By micro-reversibility, there is a one-to-one correspon- 
dence between each trajectory in the first situation with 
the time reversed trajectory in the second situation. Fur- 
thermore, since Hamiltonian dynamics preserves phase 
volume, the numerators in the r.h.s. of Eqs. (|Bip and 
(IB2I) are identical. We conclude that 



P (NitPiiUi) (AN,Ap,AU) 
P (N '- P 'u!)(-&N,Ap,-AU) 



n 



n 



exp( 



(Ni, Pi ,Ui) 

,AS\ 



kr 



(B3) 



where we used the fact that the entropy of a state is 
Boltzmann's constant time the logarithm of the phase 
space volume fi of that state. AS is thus the en- 
tropy difference between states with and without the 
primes. We now note that inverting the momenta of 
the gases does not change the statistics of particle and 
energy transport, but will obviously change the sign of 
the momentum transfer, P(n' ,- P '.,u')(— AN, Ap, —AU) = 
Pcn'.,p' ,u')(~ AN, — Ap, — AU). Hence we can rewrite 
Eq. ' ([B3)) 1 as follows: 

A <? 

exp( — ) (B4) 



P(N i ,p i ,u,)(AN, Ap, AU) 
p (Ni, P >,u{)(-AN,-Ap,-AU) 



Finally, we consider the thermodynamic limit of in- 
finitely large systems with finite particle density rii = 
Ni/Vi, momentum density mVi = Pi/Ni, and energy den- 
sity Ui = Ui/Vi- We furthermore assume that the effects 
of the removal of the constraint do not scale with the 
volume, so that during the finite time t it results in non- 
extensive changes in the parameter values. Hence we 
can drop the sub-indices of P on the l.h.s. of Eq. (|B4I 
Adding the sub-index t to emphasize the duration of the 
exchange, one can thus rewrite Eq. (|B4[) as Eq. ([7]) of the 
main text. 



APPENDIX C 

Here we briefly derive the formula for the transition 
rate T A ^ B (E,p x ) using the kinetic theory of gases. We 



7 



consider the z direction to point from reservoir A to reser- 
voir B. We require TA^B(E,p x )dEdp x dt, the probabil- 
ity to observe a particle with kinetic energy in the range 
(E, E + dE) and momentum in the range (p x ,p x + dp x ) 
to cross the hole from A to B in a time interval dt. The 
z component of the position of a particle with velocity 
v must be located within a cylinder of base area a (the 
area of the hole) and height v z dt measured from the wall. 
Furthermore, it must be traveling in the +e z direction 
(towards the hole). The appropriate expression is 



T a ^b(E,p x ) = 

dv x I dv y I dv z av z n A (j) A {v,VA) 



where 4>a(v, Va) is the Maxwellian given by Eq. ([2]), and 
we have explicitly noted its dependence on v and Va- 
A similar equation holds for Ta—>b- The v x integral is 
trivial because of the second delta function, and the re- 
maining integrals can be easily performed by changing to 
polar coordinates (R,9), given by 



R z 



tan 8 = 



(C2) 



(CI) This yields the expressions given in Eq. (fTTj) . 



x (5( — E)S(mv x -p x ), 
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